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ANNUAL REPORT 

Technical Proposal 

The objective of the proposed work is to continue the work to de- 
velop, verify and incorporate the baseline two-equation turbulence mod- 
els which account for the effects of compressibility at high speeds into a 
three-dimensional Reynolds averaged Navier-Stokes (RANS) code and 
to provide documented descriptions of the models and their numerical 
procedures so that they can be implemented into the NASP CFD codes. 

Computational Fluid Dynamics (CFD) is recognized as a signif- 
icant engineering design tool in modern hypersonic projects such as 
the National Aerospace Plane (NASP). The design of the NASP ve- 
hicle is a highly complex process. One of the critical tasks involved 
in such a process is the ability of the turbulence model to predict 
hypersonic viscous/inviscid interactions, mixing problems, transition, 
chemical nonequilibria, and a range of other turbulence phenomena. 

Turbulence models are developed on the basis of insight gained 
from experimental and theoretical research. The complexity of tur- 
bulence requires that the mathematical models be guided by the flow 
physics in a rational and practical approach. Test and validation of 
new models with data of recognized quality is an essential step toward 
model acceptability. A database with a detailed experimental infor- 
mation on high speed turbulent flows is currently being completed. 
This database provides a benchmark for testing and validation of the 
compressible model modifications. Specific topics for the database are 
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shock wave/boundary layer interactions, high-speed attached boundary 
layers with pressure gradients, and supersonic shear layer mixing. 

The baseline k-u and k - e models with the proposed length-scale 
and rapid compression modifications have already been implemented 
into the 3-D Navier-Stokes base code. The turbulence models have 
been tested against experimental data of the hypersonic database col- 
lection. The validation ranges from simple flat plate flows to complex 3- 
D flows. They included theoretical correlations, Mach 5 supersonic flat 
plate turbulent flow, Mach 3 supersonic compression corner with flow 
separation, and Mach 8 hypersonic 3-D shock-wave boundary layer in- 
teraction flows generated with single fins and two oblique fins mounted 
on a flat plate surface. A major consideration throughout the research 
effort will continue to be the development of improved compressibility 
corrections to the turbulence models and the identification of models 
which are superior in their predictive capabilities. 
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Work Statement 


The following is the work statement proposed to be accomplished 
within this past research year: 

1. Incorporate k - e models into base code and evaluate performance 
against 2-D compressible turbulent flows (Mach 5 flat plate flow 
and Mach 3 compression ramp). 

2. Test and evaluate improved numerics that take into account second 
law of thermodynamics and higher-order TVD flux limiters. 

3. Study and continue development of compressibility corrections for 
3-D turbulence models (numerical simulations of 3-D shock wave / 
boundary layer interaction (SWBLI) at Moo = 8.2). 

4. Supervise and deliver information needed to implement compress- 
ible two-equation turbulence models into 3-D codes. 
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Work Accomplished 

1. The baseline turbulence models have been incorporated into the 
base RANS 3-D Reynolds Averaged Navier-Stokes code and vali- 
dated against the Mach 5 flat plate boundary layer flow and the 
2-D codes of T.J. Coakley and RG. Huang. The two equation 
turbulence models already implemented are: 

• Wilcox k-u> model (baseline model) 

• Wilcox transition k-u model 

• Jones-Launder k - e model 

• Launder-Sharma k - e model (baseline model) 

• Chien k - e model 

2. Four important model modifications have been implemented and 
subject to different database tests: 

• A length scale modification in order to decrease the heat trans- 
fer rate in flow reattachment zones (all models). 

• A rapid compression modification to increase the size of the 
flow separation zone in high speed flows (all models). 

• A rotation or vortex stretching length scale modification to de- 
crease the level of turbulent dissipation in flow recirculating 
zones (all models). 

• A turbulence transition modification to simulate flow transition 
and complex surfaces with both laminar and turbulent com- 
pressible flow zones (k — w model). 

3. Numerical simulations with and without the different combinations 
of model modifications have been done. These simulations have 
been based on the following flows: 

• Mach 14 laminar boundary layer on a 24 ° compression ramp 
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(code validation). 

• Mach 5 flat plate turbulent boundary layer (turbulence model- 
ing validation). 

• Mach 3 separated turbulent boundary layer on a 24 ° compres- 
sion ramp (turbulent flow separation validation). 

• Mach 8 shock wave / turbulent boundary layer interaction gen- 
erated with a 15 ° and a 10 ° fin mounted on a cooled flat plate 
with T wa ii/Too = 0.3 and Re 00 = 5- 10 6 (complex 3-D SWBLI with 
crossflow separation). 

• Mach 8 intersecting shock waves / turbulent boundary layer 
interaction generated with two 15 ° fins mounted on a flat plate 
(complex 3-D SWBLI with completely crossflow separation). 

4. Important results: 

• The baseline models provide excellent agreement with the ex- 
perimental data in the prediction of the pressure profiles, skin 
friction, and heat transfer rate on the cooled flat plate of the 
complex 3-D SWBLI flow at Mach 8. 

• The baseline models provide good agreement with the experi- 
mental measurements of pitot pressure in the crossflow zone of 
the 3-D SWBLI. 

• The length scale and rapid compression modifications do not 
improve the prediction of the effects of crossflow separation. 
Although, they have been shown to significantly improve the 
predictions of 2-D hypersonic separated flows. 

• The standard models provide good agreement with experimen- 
tal measurements of flow separation in the Mach 3 supersonic 
flow ( 24 ° compression ramp). However, they underpredict the 
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size of flow separation in hypersonic flows at Moo = 7 and Moo = 9 
(P.G. Huang and T.J. Coakley) and transonic flows (F . Menter). 

• The rotation modification together with the length scale and 
rapid compression modifications improve the prediction of flow 
separation, while the length scale and rapid compression mod- 
ifications overpredict the size of flow separation in the Mach 3 
supersonic flow. 

• The models do not have the capability to predict the heat trans- 
fer rate and the skin friction on the laminar fin with a turbu- 
lent boundary layer on its junction with the flat plate. All the 
models overpredict the experimental measurements. This is a 
weakness of the present models that we plan to analyze. 

• The low Reynolds number transition modification proposed by 
Dr. David Wilcox improves the agreement with experimental 
data on transition Reynolds number in the Mach 5 flat plate 
bo un dary layer. It predicts poorly the heat transfer rate on 
the fin surface (turbulent / laminar boundary layer) of the 3-D 
SWBLI flow. 

• The transition Reynolds number is very sensitive to the inflow 
turbulence levels and the free stream turbulence. The model is 
also very sensitive to freestream turbulent kinetic energy (Dr. 
Thomas Coakley). Parallel research has shown that small val- 
ues of free stream specific dissipation rate u has strong influence 
in the prediction of turbulent free shear flows and weaker influ- 
ence in the prediction of attached boundary layers (F. Menter). 
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5. Improved numerics have been tested and implemented into the 3-D 
RANS Code. 

• Different higher-order finite difference numerics and TVD flux 
limiters methods have been studied and developed for the nu- 
merical method of the turbulence equations and the Reynolds- 
Averaged Navier-Stokes equations. 

• This research has shown that there axe unlimited number of 
averaging procedures for flux-difference splitting methods that 
are consistent with Roe’s conservative property U. In general, 
the Jacobian matrix of the discrete system is composed of a 
particular solution and an homogeneous solution. Roe’s geo- 
metric averaging procedure represents one of these solutions, 
while the arithmetic averaging procedure incorporated in the 
RANS3D code represents another solution. All these methods 
are consistent in their limit with the Jacobian matrix of the 
continuos system. 

• The second-order flux limiters and the flux-difference splitting 
based on arithmetic averaging have been improved and tested 
in the simulation of the M = 14.1 laminar hypersonic 24° ramp 
and the turbulent flow at M = 5 on the cooled flat plate. The 
results are consistent with previous simulations and the running 
time improved 18%. 

• The code has been improved in user friendliness with a general 
input subroutine that reads the input data (parameters of tur- 
bulence models, different implicit boundary approximations in 
each direction, etc..), instead of using an specific code for each 
case (such as the CNS code). 


7 



• The 3-D code has been delivered to Dr. Florian Menter (Model- 
ing and Experimental Validation Branch) in order to study the 
ability of the turbulence models and modifications to predict 
transonic flows. It has been requested by Dr. Ethiraj Venkat- 
apathy (Aerothermodynamics Branch) in order to study shock 
waves / turbulent boundary layer interactions in inlets. It has 
been requested by McDonnell Douglas Space Systems Company 
to study high-speed turbulence aero-optics. 

• Plotting programs for the 3-D simulations based on PV-WAVE 
coupled with FRAME software have been developed. These 
programs reduce and plot different combinations of mean vari- 
ables and turbulence statistics. 

• The numerical simulations have been rim with the Cray-YMP 
and the Cray-2 NAS computers at NASA Ames Research Cen- 
ter. 

6. A paper 1 has been written and presented with different results of the 
3-D Navier-Stokes simulations in the AIAA Fourth International 
Aerospace Planes Conference held on December 1992 in Orlando, 
Florida. The paper is titled “Two-Equation Turbulence Modeling 
for 3-D Hypersonic Flows,” and co-authored with T.J. Coakley and 
J.G. Marvin. The organizer of the Meeting session was Dr. John 
Hicks, NASP Chief Engineer. 
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Future Work 


1. Develop a theoretical study of the improved compressibility cor- 
rections to the turbulence models in the boundary layer equations 
to identify their superior predictive capabilities. Present models 
do not support the ‘production equals dissipation’ balance in high 
speed flows. The model equations need to be analyzed by wall 
zone layers with appropriate boundary interface of the nonlinear 
equations. 

2. Test and evaluate refinements of the baseline two-equation turbu- 
lence models implemented in the previous phase. Research focused 
in compressibility corrections for flow separation, transition, and 
low Reynolds number coefficients in wall sublayers. 

Incorporate the following low Reynolds number modifications 
into base code 2-4 : 

• Huang-Coakley k - e model. 

• So-Zhang-Speziale k - e model. 

• Coakley q - u model. 

• Menter k-u/k - e model. 

• Rodi k - e with one equation model in sublayer zone. 

Validate the turbulence models against experimental data of 
the following flows 5 : 

• Mach 5 flat plate flow. 

• Mach 3 compression corner flow. 

• Mach 9 compression corner flow. 

• Mach 7 cylinder-flare flow. 

3. Study and development of baseline modifications for 3-D turbulence 
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flows: 

3a. Simulate 3-D SWBLI (1 fin), Moo = 8.2, of M.I. Kussoy and K.C. 
Horstman 6 . 

3b. Simulate 3-D intersecting SWBLI (2 fins), M « = 8.3, of M.I. 
Kussoy and K.C. Horstman 7 . 

3c. Simulate 3-D intersecting SWBLI (2 fins), Moo = 4, of T.J. Gar- 
rison and G.S. Settles 8 . 

3d. Simulate 3-D intersecting SWBLI (2 oblique fins), Moo = 8.3, of 
Bahrami, P. et al 9 . 

3e. Simulate 3-D intersecting SWBLI (2 fins with top wall), Moo = 
8.3, of M.I. Kussoy and C.C. Horstman 10 . 

4. Documentation. Supervise and deliver information needed to use 
the 3-D RANS code, to implement the baseline turbulence models 
into the NASP 3-D codes, and to documented descriptions of the 
models and their numerical procedures so that they can be imple- 
mented into the NASP CFD codes. 

5 3-D Navier-Stokes code. Improve code with self adjusting parame- 
ters to run multiple 3-D grids with multiple segments and multiple 
boxes, main effort to simulate flows with complex physics to im- 
prove understanding of the turbulence effects. 
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Adiabatic Flat Plate 


® Atoo — 5* 

• Ree = 10 4 . 

• Turbulence models: 

• k - e Launder-Sharma. 

• k - to Wilcox. 

• k — u) baseline (includes algebraic length scale limit and rapid 
compression modifications). 

• k - x with rotation modification. 
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VELOCITY PROFILE, M=5 
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Supersonic Compression Corner 


• Settles 4 24° compression corner, M = 2.84. 

• Turbulence models: 

• k -oj Wilcox. 

• k - u> baseline (includes algebraic length scale limit and rapid 
compression modifications). 

• k - oj wlr-mod (includes baseline and rotation modifications). 

• k-oj 1-mod (includes algebraic length scale limit modification). 

• k - oj r-mod (includes rapid compression modification). 
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ABSTRACT 

An investigation to verify, incorporate and de- 
velop two-equation turbulence models for three- 
dimensional high speed flows is presented. The cur- 
rent design effort of hypersonic vehicles has led to an 
intensive study of turbulence models for compressible 
hypersonic flows. This research complements an ex- 
tensive review of experimental data and the current 
development of 2-D turbulence models. The review of 
experimental data on 2-D and 3-D flows includes com- 
plex hypersonic flows with pressure profiles, skin fric- 
tion, wall heat transfer, and turbulence statistics data. 
In a parallel effort, turbulence models for high speed 
flows have been tested against flat plate boundary lay- 
ers, and are being tested against the 2-D database. 
In the present paper, we present the results of 3-D 
Navier-Stokes numerical simulations with an improved 
k — two-equation turbulence model against experi- 
mental data and empirical correlations of an adiabatic 
flat plate boundary layer, a cold wall flat plate bound- 
ary l aver, and a 3-D database flow, the interaction of 
an oblique shock wave and a thick turbulent bound- 
ary layer (SWBLI) with a free stream Mach number 
Moo — 8.18 and Reynolds number Re a Q = 5 • 10 6 . 

INTRODUCTION 

The cost-effective engineering design of new vehi- 
cles encountering subsonic, transonic, supersonic, and 
hypersonic speeds requires an advanced and efficient 
computational fluid dynamics (CFD) technology for 
the accurate aerodynamic prediction of the complex 3- 
D flow fields. The CFD development as a design tool 
of complex vehicles with complex flow fields is based 
on advances in the understanding of the physics of the 
flows, experimental observations of turbulent flows, 
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advances in numerical methods, and availability of 
faster numerical processors with large memory capac- 
ity. This method is needed to accurately predict the 
critical aerodynamic loads and wall heating rates in 
the design of hypersonic vehicles. The combination 
of an efficient three dimensional CFD shock capturing 
method with an accurate turbulence model is here re- 
quired to account for the effects of the interaction of a 
planar shock wave which sweeps across the turbulent 
boundary layer. If the fin angle is large, the wall pres- 
sure rise and the boundary layer cause the formation of 
a strong crossflow separation vortex below the inviscid 
shock wave and inside the turbulent boundary layer. 
Part of the flow over the separation vortex reattaches 
behind the vortex impinging on the plate and leading 
to larger increments of wall pressure, skin friction, and 
heat transfer. 

The interaction of shock waves with turbulent 
boundary layers causes considerable efficiency losses 
and represents a major design problem of vehicles fly- 
ing at hypersonic speeds. This kind of interactions 
determines the performance of inlets for airbreath- 
ing propulsion at hypersonic speeds. Extensive ex- 
perimental studies have been conducted for the 3-D 
sharp fin/ flat plate geometry 1 ” 7 . Shock wave/ bound- 
ary layer interactions is the first specific topic of the 
hypersonic database collection and assessment effort 
of the critical review sponsored by the NASP pro- 
gram through NASA Ames Research Center 1 . This 
database collection ha s been generated to provide a 
common base for testing hypersonic turbulence mod- 
els. This collection was based on a review of a survey 
of thousands of references which was narrowed to 105 
distinct experimental studies. The experiments were 
subject to a database assessment and a few of them 
were reported in Category I 2 (accepted) or Category 
II 3 ’ 4 (limited acceptance). The data include wall pres- 
sure, skin friction, wall heat transfer, and flow field 
surveys at supersonic/ hypersonic free stream Mach 
numbers 3 through 6, with an inviscid pressure ratio 
across the planar shock varying from 1.35 to 5.1, re- 
spectively. At this time, the database is under further 
study in order to include new hypersonic experimental 
data. 

Recently, an experimental study at Mach 8.18 was 
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conducted by Kussoy and Horstman 6-8 . This study 
includes measurements of pressure profiles, skin fric- 
tion, and surface heat transfer with 5°, 7.5°, 10°, 12.5°, 
and 15° fin angles. The experimental data also in- 
clude flow field yaw angles and pitot pressures with 
10° and 15° fin angles, corresponding to an inviscid 
pressure ratio of 5.4 and 9.7, respectively. Compu- 
tational studies by Knight, Horstman, and Monson 8 
include 3-D Reynolds averaged Navier-Stokes simula- 
tions with two k — € turbulence models and a “conical” 
simulation with the Baldwin-Lomax mixing length tur- 
bulence model. Their 3-D Navier-Stokes simulation 
shows agreement with the experimental measurements 
of surface pressure, skin friction, flow direction, and 
wall heat transfer on the plate; it failed to predict the 
initial heat transfer peak and pressure valley observed 
below the inviscid shock wave; and the peaks of the 
heat transfer rate distributions are overpredicted by 
36% and 48% in the 10° and 15° fin angle, respec- 
tively. The heat transfer rate on the fin section is over- 
predicted by a factor of 2 approximately, “due to the 
assumption of fully turbulent flow from the fin lead- 
ing edge.” The simulation with a laminar boundary 
layer on the fin overpredicts the heat transfer data 
inside the boundary layer thickness and agrees with 
the data outside of the boundary layer edge. Al- 
though the experimental boundary layer is laminar 
on the fin above the flat plate boundary layer, the 
simulations with turbulent and laminar flow on the 
fin show similar results of pressure, skin friction, and 
heat transfer on the flat plate. On the other hand, 
the “conical” simulations show poor agreement with 
experimental measurements. The wall pressure in the 
fin/ plate intersection zone is overpredicted by at least 
25%, while the peak heat transfer rate is overpredicted 
by at least 200%. These differences were ” attributed 
to the limitations of the Baldwin-Lomax turbulence 
model and/or the assumption of conical flow.” The 
present simulations shown in this paper indicate that 
their differences are mostly due to their “conical” flow 
assumption. 

Several experimental investigators 9- 11 have been 
advancing the assumption of quasi- conical symmetry 
based on the empirical and numerical observations 
that the interaction vortex grows with distance follow- 
ing the direction of the inviscid oblique plane shock 
wave. This approximation is supported by the wall 
pressure data at lower Mach numbers in zones far 
from the fin leading edge and far from the intersec- 
tion of the flat plate and the fin 11 . It is not sup- 
ported by mean surface data of skin friction and heat 
transfer rate 5,12 . The main purpose of using conical 
coordinates in numerical simulations is to reduce the 
computer time needed to solve these flows and to in- 
crease the grid resolution. In general, the quasi-conical 
simulation® required about 1.2 hours of Cray-YMP 
with a 98x88 grid and about 3.6 hours of Cray-YMP 


with a 195x191 fine grid, while the 3-D Navier-Stokes 
simulations 8,13 with the hybrid explicit-implicit algo- 
rithm of MacCormack 14 required about 25 hours of 
Cray-YMP with a 64x40x60 grid. 

Parallel to the experimental research effort, tur- 
bulence models for high speed flows are currently being 
developed and tested for the design of hypersonic ve- 
hicles. These models include different versions of the 
4 — e and 4— w two equation turbulence models 15 . The 
numerical results show good agreement within 10% of 
skin friction measurements and empirical correlations 
in adiabatic and cooled- wall flat plate boundary lay- 
ers (0.2 < T w /T a tv <1.0 and < 10). Their heat 
transfer predictions show a Reynolds analogy factor 
of 1.17 to 1.19. This kind of agreement is also ob- 
served when the proposed modeling compressibility 
corrections for complex flows are incorporated into the 
models. Several corrections have been proposed to re- 
duce the overprediction of mean heat transfer rates 
in reattachment zones of separating boundary layers, 
to increase the size of flow separation and the corre- 
sponding surface pressure distributions, and to reduce 
the spreading rate of turbulent shear layers, to im- 
prove the prediction of flow separation. These tur- 
bulence model corrections showed significant improve- 
ments in the wall pressure and mean heat transfer 
rate predictions of 2-D shock wave/ boundary layer 
interactions 16 . Numerical simulations are beginning to 
provide significant results in engineering design; in par- 
ticular, efficient computational fluid dynamics (CFD) 
methods are beginning to serve as major complement 
to wind tunnel experimental research due to the latest 
advances in CFD and available computer power. 

The objectives of this paper are to evaluate a 
baseline 4 — model in a strong shock wave/ bound- 
ary layer interaction, to examine the capability of the 
simulation to represent the complex 3-D structure of 
the flowfield, and to demonstrate the capability of a 
new efficient 3-D computational Navier-Stokes method 
with a two-equation turbulence model. 

3-D NUMERICAL METHOD 

The numerical method 17 is based on a 3-D im- 
plicit upwind flux-difference split finite difference al- 
gorithm of the mean compressible Navier-Stokes equa- 
tions using mass-averaged variables. The conserva- 
tion laws of mass, momentum, and energy for the 3-D 
Navier-Stokes equations are expressed in generalized 
curvilinear coordinates and in compressed vector no- 
tation as 

(JU) it +7t~F jt i = 0 ( 1 ) 

where repeated sub-indices in any term imply sum- 
mation over the index range, sub-index t following a 
comma implies partial differentiation with respect to 
time, sub-index j following a comma implies partial 
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differentiation with respect to the spatial Cartesian co- 
ordinate direction Xj , sub-index i following a comma 
implies partial differentiation with respect to the spa- 
tial curvilinear coordinate direction 4» , and the overbar 
denotes spatial averaging. The vector U represents the 
conservative dependent variables, 

U = (p,pui,pti 2 ,/)ti3,e) (2) 

Fj are the flux vectors in the respective Cartesian co- 
ordinates direction Xj , and J is the Jacobian of the 
coordinate transformation. The metric coefficients are 
evaluated at the interface between the grid points of 
the spatial difference. 

All flux differences are treated implicitly in order 
to increase stability and to be able to use large incre- 
ments of time or CFL numbers. The numerical scheme 
for the viscous fluxes is second-order central difference, 
while the numerical scheme for the inviscid fluxes is 
a higher-order TVD upwind flux-difference splitting. 
Different tests and experiments in subsonic, transonic, 
and supersonic flows have shown that the conservation 
properties are maintained once convergence is achieved 
with a second- or higher-order scheme 18,19 . 

The differentiation of the inviscid fluxes is similar 
to Roe’s flux difference splitting method, however it 
is based on simple arithmetic averaging of primitive 
variables. In order to preserve the conservative prop- 
erty of the fluxes, each flux difference term ■ of the 
conservation law equations is defined as 

8F = F- +1 - Fj = A SU = SAS- 1 SU (3) 

where the subindex j is omitted in equation 3 for sim- 
plicity. A is a diagonal matrix whose coefficients are 
the eigenvalues of the matrix A, and 5 is the column- 
eigenvector matrix of A. The inviscid first-order flux 
differences terms are split according to the sign of 
the eigenvalues and unconditionally stable Euler im- 
plicit methods are constructed by backward upwind 
differencing flux differences with positive eigenvalues 
and forward differencing flux differences with negative 
eigenvalues. 

Higher-order spatial TVD flux differences in the 
right-hand-side of the inviscid terms of the conserva- 
tion law equations are defined by using the n minmo<T 
limiter of Osher and Chakravarthy 20 applied on the 
complete split flux difference instead of the charac- 
teristic variable difference. It has the capability to 
represent various higher-order differences: first-order 
upwind, second-order upwind, third-order upwind bi- 
ased, second-order Fromm scheme, and other combi- 
nations of upwind and central differences. 

The solution procedure is based on an implicit 
” method of planes” symmetric Gauss-Seidell relax- 
ation scheme. The data is conveniently stored on suc- 
cessive planes along the streamwise coordinate, and 


the system of equations is solved in each successive 
plane along the forward direction, first, and along 
the backward direction, afterwards In each plane, 
the solution is obtained by using a two level pseudo 
time dependent relaxation procedure based on a di- 
agonally dominant approximate factorization DDADI. 
The space marching alternating directional sweeps in 
the streamwise direction are von Neumann uncondi- 
tionally stable for zones of subsonic and streamwise 
separated and reversed flows as well as supersonic flow. 
The space marching method results in improved prop- 
agation of nonlinear effects to accelerate convergence 
to steady state, much as do the more restrictive PNS 
techniques 18,19 . 

The diagonal dominant approximate factorization 
of the left-hand-side of the conservation law equations 
including the implicit viscous terms leads to the fol- 
lowing block tridiagonal equation sequence for the 
plane relaxation method 

(-A+, D, A-J8U' = - RHS "’ n+1 (4a) 

(~A+ 3 , D, Af 3 )6U --D8U’ (46) 

The diagonally dominant matrix D involves the first- 
order split Jacobian matrices and the Jacobian matri- 
ces of the viscous terms of each coordinate direction 
and the solution is updated from time step n to time 
step n + 1. Observe that the RHS of equation 4a has 
an exponent n, n-h 1 because some terms in the stream- 
wise direction are already updated at time step n+1 
due to the plane relaxation procedure. 

A Newton-Raphson acceleration procedure is ob- 
tained by solving each plane twice or more times in 
each relaxation sequence. This procedure produces 
significant improvement in the propagation of nonlin- 
ear effects due to the nonlinearity of the Jacobian ma- 
trix. 

BOUNDARY CONDITIONS 

Mathematically well posed implicit boundary con- 
ditions are imposed in every grid boundary point. 
These conditions are mathematically well posed and 
based on the sign of the eigenvalue of the correspond- 
ing component of the flux-difference split vector. The 
physical boundary approximations imposed on the cal- 
culations are frozen inflow, no slip velocity and adia- 
batic or cold temperature on the wall surface, symme- 
try plane, and flow through in free boundaries. The 
inflow profiles are defined by the free stream values on 
the flat plate calculations, and on the flat plate profiles 
with matching 8* experimental values on the fin/flat 
plate calculations. At the flow through and outflow 
boundary planes, if the flow is subsonic inflow then 
no change in entropy, enthalpy and tangential veloc- 
ity are imposed, and if the flow is subsonic outflow 
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then no pressure variation is imposed in the differenc- 
ing toward the boundary. Free stream static pressure 
is effectively maintained, and outflow boundary data 
is properly extrapolated by the algorithm. The sym- 
metry plane condition is applied to the grid points on 
the plane upstream of the fin leading edge. 

TURBULENCE MODELS 

Two kind of turbulence models have been used in 
these simulations. The first kind of model is the k-u> 
two-equation turbulent eddy viscosity model 15 . The 
eddy viscosity is defined as 

pr — Cp p k/uf (5) 

where k and w are the additional field variables. Ob- 
serve that u> is defined a s e/k whereas the traditional 
formulation is t/c^k. The k — u> field equations are 
written as 

(Jpk) ti + J£i t j ( pkuj )■ = J(c ki P - Ckipuk) 

- (JZiimiP + ( 6 ) 

{Jpw),t + j£i ,j (putuj)j = J{cu>\Pu/k - C w 2 p^ 2 ) 

— d* CO 

where P the production of turbulent kinetic energy is 

P - +uj ti )uij - 3) - pkuij 2/3 (8) 

and the model constant or parameters are 

an = 1.0 c wl = 0.555 & k — 2.0 = 0.09 

c M =1.0 c w2 = 0.833 <7* = 2.0 (9) 

This model does not need explicit damping terms in 

the model parameters due to the behavior of u> in the 
wall layer, although Wilcox 21 has proposed a low tur- 
bulent Reynolds number dependence of the model pa- 
rameters (c^, Cfc 2 , and c w i) in order to describe tran- 
sition in simple 2-D flows. The boundary conditions 
at the wall are k = 0 and u) = $p/(c u/ 2 py\) where 
yi is the wall distance from the first grid point off 
the wall. No difference is found if w is redefined after 
each iteration by the minimum value between u> itself 
and 6/i/(c u ,2/>y 2 ) within the laminar sublayer, however, 
there is an improvement in convergence speed. 

The modeling correction for reduction in the pre- 
diction of heat transfer rate in the reattachment zone 
is based on an algebraic upper limit of the turbulent 
length scale near the wall 15 . This limit is defined by 
the Prandtl length scale, «y, and it is used to redefine 
the specific dissipation rate u = Vk/l by the following 
relationship, 

/ = min(c; 3 / 4 «y, y/k/u) and /c = 0.41 (10) 


This modification is mainly significant near reattach- 
ment and has little effect in other zones 15,16,22 . The 
modeling correction to increase the separation bub- 
ble size is based on the assumption that the prod- 
uct pi remains constant under rapid compression and 
it adds a term proportional to the dilatational term, 
—2pkui } i/3, -72/3(2 - c Wl )puni i>t to the right-hand- 
side of the (j eq. 7. In other words, the net contribu- 
tion of the dilatational term, including the last produc- 
tion term of Eq. 8, is increased from — /2/3c u , 1 /xju t| i 
to -J4/Zpuui t i. The k-u> model including the length 
scale and the rapid compression modifications is called 
the baseline k — w model. 

The second kind of model is an algebraic turbu- 
lent eddy viscosity model based on the Baldwin and 
Lomax 23 formulation. This model is a two-layer mix- 
ing length formulation for the eddy viscosity similar 
to the mixing length model of Cebeci 18 . The inner 
layer uses the Prandtl- Van Driest formula with a lower 
bound value of 1/14 of the maximum viscosity value in 
free shear zones, while the outer layer uses a function 
of the maximum value of the mean vorticity and total 
velocity difference. 

DESCRIPTION OF THE SIMULATIONS 

The computational grids have been generated us- 
ing simple algebraic cubic interpolation and clustering 
to the wall boundaries. The number of grid points 
in these simulations are 3x61x61 for the Mach 5 flat 
plate flow (see Fig. la), 91x61x3 for the Mach 8.18 
flat plate flow (see Fig. 2a), and 61x41x61 for the 
Mach 8.18 3-D oblique shock wave - turbulent bound- 
ary layer interaction flow (see Fig. 3). The Mach 5 
flat plate grid has 3 parallel planes to the mean flow, 
and the solution procedure applied was the diagonally 
dominant approximate factorization in the midplane. 
The value of y + for the first grid point off the wall 
varied from 0.2 to 0.35, and about 45 grid points were 
located inside the boundary layer. The Mach 8.18 flat 
plate grid has 91 planes along the mean flow, and the 
solution procedure was the plane relaxation with four 
inner Newton-Raphson iterations. The value of y + for 
the first grid point off the wall was 0.02, and about 
40 grid points were located inside the boundary layer. 
This grid was clustered to the wall and to the lead- 
ing edge in order to provide the initial profiles for the 
shock wave - boundary layer interaction (SWBLI). The 
61x41x61 grid is similar in topology to the 64x40x64 
grid used by Knight et al 8 in a previous 3-D SWBLI 
computation. Each longitudinal grid line along the flat 
plate has 61 grid points from the inflow through the 
outflow plane; each grid line normal to the flat plate 
has 41 grid points from the plate through the outer- 
flow plane; and each grid line normal to the mean flow 
direction and the flat plate has 61 grid points from 
the fin through the free flow lateral plane. A coarser 
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31x21x31 grid which uses every other grid point of the 
finer 61x41x61 grid was used to achieve faster con- 
vergence time, to examine different model corrections, 
and to examine grid effects in this simulation. The 
solutions of the coarse grid were used to provide the 
initial flowfield of the fine grid. Figure 3 shows a 3-D 
view of the grid points located on the flat plate and 
the fin surface. The minimum grid point spacing nor- 
mal to the walls has been fixed to 4 ■ 10“ 6 m; in wall 
coordinates, the grid point spacing at the inflow sec- 
tion of the plate is y + = 0.08. the solution procedure 
used the plane relaxation in the streamwise direction 
and the diagonally dominant approximate factoriza- 
tion along the curvilinear coordinates on the plane. 

The 3-D Navier-Stokes simulations of the SWBLI 
with the coarser grid required between 30 to 50 min- 
utes, while the simulations with the finer grid required 
between 3 to 6 hours of CRAY-YMP cpu time. The 
adiabatic and cooled-wall flat plate simulations re- 
quired only a few minutes of cpu time. 

All the simulations assume perfect air with con- 
stant ratio of specific heat y = 1.4, molecular Prandtl 
number 0.73, and turbulent Prandtl number 0.9. The 
laminar viscosity is specified by Sutherland’s law 24 . 

ADIABATIC FLAT PLATE FLOW, M = 5 

The main objectives of the adiabatic flat plate cal- 
culations were to verify the 3-D numerical simulations 
and the turbulence models. Following the detailed 
comparison of Coakley and Huang 15 in the develop- 
ment and testing of two-equation models in 2-D flows, 
a comparison with empirical correlations of experimen- 
tal data and other numerical simulation is presented. 
The mean flow conditions are Moo = 5, T w /T aw — 1, 
and Reoo = 2.2- 10 7 per meter. This Reynolds number 
avoid low Reynolds number effects 4 . Figure lb shows 
the comparison of the van Driest II empirical correla- 
tion, the skin friction results using the van Driest com- 
pressible transformation 25 shows excellent agreement 
with the Karman-Schoenherr correlation curve 26 as a 
function of the momentum thickness Reynolds num- 
ber. Figures lc and Id show the comparisons of the 
velocity profile with the compressible law of the wall 25 
in wall coordinates, and the 2-D prediction of Coak- 
ley et al 15 , respectively. The results show very good 
agreement in the prediction of the law of the wall and 
velocity profile obtained with a different 2-D code 15 . 
As it is expected, there is little difference between the 
results obtained with the k — lj model with and with- 
out the proposed modeling corrections described in the 
turbulence models section. 

COLD FLAT PLATE FLOW, M^ = 8.18 

The cold flat plate flow calculations were done to 


simulate the experimental results of Kussoy et al 4 , in 
order to provide the initial profiles of the SWBLI sim- 
ulations described below, and to verify the 3-D nu- 
merical simulations. The mean flow conditions are 
Moo = 8.18, T w /Tx oo = 0.27, T w /T aw — 0.3, and 
Rtoo — 5 • 10 6 per meter. Following the detailed 
comparison of Coakley and Huang 4 in the develop- 
ment and testing of two-equation models a comparison 
with empirical correlations of experimental data and 
other numerical simulation is presented at Figure 2b 
shows the comparison of the van Driest II empirical 
correlation, the skin friction results using the van Dri- 
est compressible transformation 25 with the Karman- 
Schoenherr correlation curve 26 . The results obtained 
with the k — u; with and without the length-scale and 
the rapid compression modifications are about 10% 
higher than those results predicted with the empir- 
ical correlation, while the results obtained with the 
Baldwin-Lomax algebraic model tend to agree with 
the empirical correlation at higher momentum thick- 
ness Reynolds numbers. The differences in the results 
may be smaller at higher Reynolds numbers. In order 
to test the predictions at higher Reynolds numbers, 
the flat plate was extended from the original length of 
2.2 meters upto a length of 3 meters in the numeri- 
cal simulations, However, the transformed momentum 
thickness Reynolds number was still less than 3000. 
Figure 2c shows the velocity profile comparison with 
the compressible law of the wall 25 in wall coordinates. 
These predictions are also similar to the 2-D predic- 
tion of Coakley et al 15 and they underpredict the law 
of the wall. Figure 2d show the full cold wall velocity 
profiles obtained with the different models. The main 
differences in the profiles are observed in the wake re- 
gion of the velocity distributions. Figure 2e show the 
typical growth of the Mach contour lines inside the 
turbulent boundary layer obtained with the modified 
k — u model, the growth rate of each contour line is al- 
most linear and uniform. Figure 2f shows the pressure 
contour levels normalized with the free stream pres- 
sure and the main feature is the pressure rise caused 
by the lip shock. 

OBLIQUE SHOCK WAVE - TURBULENT 

BOUNDARY LAYER INTERACTION 

The experiment was conducted with free stream 
Mach number Moo = 8.18, temperature = 81 K, 
and pressure poo = 430iV/m 2 . The flat plate dimen- 
sions were 2.2 m length and 0.76 m width, while the 
fin dimensions were 0.3 m length and 0.2 m height. 
The fin was mounted 1.76 m downstream of the flat 
plate leading edge. Both the fin and flat plate were 
cooled with a temperature of 300 ° K which is equiva- 
lent to T w = 0.27 T too = 0.3 T aw . The free stream Mach 
number is 8.18 and the free stream Reynolds number 
is 5 • 10 6 per meter. 
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The inflow profiles for the 3-D SWBLI simulations 
were obtained from the flat plate profiles that matched 
the experimental value of the compressible displace- 
ment thickness, 5 * =: 0.0159m. The experiment was 
done on a flat plate with -2° angle of attack in order 
to attain a uniform 2-D inflow before the fin interaction 
zone. 3-D Navier-Stokes simulations of the flat plate 
flow were run with and without angle of attack. As 
it should be expected, differences are mainly obtained 
in the magnitude of the primitive variables inside the 
shock layer. Since most of the experimental results 
are nondimensionalyzed with shock layer values, the 
inflow boundary layer profiles were generated by run- 
ning a Navier-Stokes simulation of the flat plate flow 
from the leading edge with no angle of attack. The 
following Figures 4 through 6 show the velocity, tem- 
perature, and density profiles, respectively, of the 3-D 
flow obtained with the k — u turbulence model includ- 
ing the modeling corrections. The numerical results 
are shown in solid lines while the experimental results 
are shown with diamond symbols. There are small dif- 
ferences in the boundary layer profiles caused by the 
difference in the angle of attack of the flat plate. 

Results for a = 10° and 15° 

The 3-D Navier-Stokes SWBLI simulations were 
done with a 10° and a 15° fin angle. Both geome- 
tries were simulated with both the coarser 31x21x31 
and the finer 61x41x61 grids. In the experiment, the 
boundary layer on the flat plate was thick and turbu- 
lent, while the boundary layer on the fin surface was 
laminar above the flat plate boundary layer edge. The 
estimated transition location of the fin boundary layer 
was downstream of the test section of the experimen- 
tal data 8 . The fin surface was simulated as a laminar 
boundary layer and also as a turbulent boundary layer, 
the results show only significant differences in the pre- 
diction of the heat transfer rate on the fin surface; this 
observation is also supported by previous simulations 8 . 
Results obtained with four different turbulence models 
are reported in this paper: the baseline k—uj , the stan- 
dard k — u) t the transitional k — lj model 21 , and the 
two-layer algebraic Baldwin-Lomax turbulence mod- 
els. Both geometries were simulated with the baseline 
k — uj and the standard k — u models. In addition, 
the flat plate with the 10° fin angle was also simu- 
lated with the Baldwin-Lomax turbulence model and 
the transitional k — u model 21 . The transitional k — uj 
model was used to test its predictive capabilities of the 
heat transfer rate data on the fin surface. 

In the following figures, the symbols show the ex- 
perimental data, the vertical lines on the symbols show 
the experimental uncertainty of the data, and the dif- 
ferent lines show the numerical results, x is the dis- 
tance on the flat plate from the fin leading edge along 
the main freestream flow, s = x/cosa is the distance 
along the fin/plate junction from the fin leading edge, 
the abscissa z — xtana represents the distance from 
the fin along the crossed section of the flat plate, and 
y represents the vertical distance along the surface of 


the fin. Figures 5 and 6 show the comparison of the 
results obtained for the 10° and the 15° fin angle, re- 
spectively. 

• The Plate Surface Pressure 

Figures 5a and 6a show the comparison of the sur- 
face pressure profile against the experimental data on 
the flat plate surface atx = 0.18I9masa function of 
distance from the fin, for the 10° and 15° fin angles, 
respectively. The ordinate P/ Poo shows the nondi- 
mensional pressure values which are nondimensional- 
ized with the freestream pressure. All the pressure 
profiles show the rise of the pressure in the bound- 
ary of the separation zone of the crossflow vortex, a 
plateau in the separation zone, and a peak pressure 
in the reattachment zone behind the vortex and very 
near the fin. The agreement between the predictions 
of the baseline k — u, the standard t - w, and the 
Baldwin-Lomax models and the experimental data is 
quite good, especially in the 10° fin angle case within 
the 10% experimental uncertainty of the data. There 
is only an underprediction of the peak reattachment 
pressure in the 15° fin angle case which it may be 
caused by lack of grid resolution and/or turbulence 
modeling. Although, some improvements in the agree- 
ment with the experimental data were obtained with 
the 61x41x61 simulations, in particular, the low pres- 
sure values shown between the pressure plateau and 
reattachment peak. The results of the two-equation 
models are closer to the peak pressure data. The re- 
sults with the baseline model are only slightly closer to 
the experimental data than those shown by the stan- 
dard k— u) turbulence model. Furthermore, the results 
with the finer grid follow the peaks and valley of the 
pressure profile more closely, which may indicate that 
grid resolution may give an even better agreement. 
The pressure profiles on the flat plate are similar to 
the predictions of Knight et al 4 using a k — e model; 
however, their “conical” solutions with the Baldwin- 
Lomax model are quite inferior and overestimate the 
peak pressure by 26% and 53% in the 10° and 15° 
fin angle, respectively. These differences should be at- 
tributed to their assumption of “conical” flow. 

• The Fin Surface Pressure 

Figures 5b and 6b show the comparison of the 
surface pressure profile against the experimental data 
obtained at $ = 0.1834 m on the fin as a function 
of the vertical distance y from the flat plate along 
the fin surface, for the 10° and the 15° fin angle, re- 
spectively. The ordinate P/Poo shows the nondimen- 
sional pressure values which are nondimensionalized 
with the freestream pressure. All the numerical re- 
sults obtained with the different turbulence models 
and grids show excellent agreement with the experi- 
mental data, in particular the baseline k — u model 
with the finer grid. The pressure profiles show a peak 
in the fin/plate junction due to the reattachment jet 
stream, a decrease inside the turbulent boundary layer 
which may be a secondary effect of the jet stream be- 
hind the main crossflow vortex, and a rise upto a high 
pressure plateau due to the pressure increase caused 
by the oblique shock wave. The pressure profile pre- 
dictions on the fin surface are a little better than those 
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presented by Knight et al 4 using a k - e model, and 
much better than their “conical” solutions with the 
Baldwin- Lomax model which are quite inferior. 

• The Plate Skin Friction 

Figures 5c and 6c show the comparison of the 
surface skin friction profile against he experimental 
data on the flat plate at i = 0.155 m as a function 
of distance from the fin, for the 10° and the 15° fin 
angle, respectively. All the 3-D numerical simulations 
show good agreement with the experimental data. The 
peak skin friction in the reattachment zone is predicted 
within the uncertainty of the data. The predictions 
are similar to those presented by Knight et al 4 using 
a k - € model, and much better than their “conical” 
solutions with the Baldwin-Lomax model which over- 
estimate the peak value by 225% and 129% in the 10° 
and 15° fin angle cases, respectively. All the numeri- 
cal predictions show a decrease in the skin friction very 
near the fin/plate corner zone. This behavior which is 
not evident from the experimental data is caused by 
a very small zone of cooled dead flow in the 90° cor- 
ner walls. The predictions do not show any evidence 
of a second smaller local peak as indicated by the few 
available data. Further improvements in turbulence 
modeling is required to represent these effects of the 
vortical structure. The predictions of the baseline and 
the standard k — u> model are very similar. 

♦ The Plate Heat Transfer Rate 

Figures 5d and 6d show the comparison of the 
surface heat transfer rate profile against the experi- 
mental data on the flat plate at x = 0.1645 m as a 
function of distance from the fin, for the 10° and the 
15° fin angle, respectively. The ordinate Q/Q 0 o shows 
the nondimensional heat transfer rate values which are 
nondimensionalized with the freestream inflow value 
upstream of the shock interaction. The experimen- 
tal profiles show a rise along the separation boundary 
of the crossflow vortex, a plateau in the separation 
zone, and a peak value in the reattachment zone be- 
hind the crossflow vortex and near the fin. For the 15° 
fin angle, the experimental profile shows a decrease 
of the heat transfer rate between the plateau and the 
reattachment peak. The agreement between the pre- 
dictions of the baseline k — u/, the standard k — <*/, 
and the Baldwin-Lomax models and the experimental 
data is quite good. Although the two-equation mod- 
els tend to underpredict the heat transfer rate under 
the separation zone, while the Baldwin-Lomax predic- 
tion is within the experimental uncertainty in the 10° 
fin angle case. The results of the simulations with 
the baseline and the standard k — w models show rea- 
sonable agreement with the data in the 15° fin angle 
case. There is an underprediction of the peak heat 
transfer rate in the separation zone, and the decrease 
of the heat transfer rate between the plateau and the 
reattachment peak is not predicted. The peak heat 
transfer rate is predicted within the 10% experimen- 
tal uncertainty in the 10° fin angle case, and slightly 
overpredicted by 14% (4% above the uncertainty of 
the data) in the 15° fin angle case. The predictions 
are also similar to those presented by Knight et al 4 
using a k — € model, however, their overpredictions of 


thepeak heat transfer rate were much higher, 36% and 
48%, respectively. Their “conical” solutions with the 
Bald will- Lomax model are also quite inferior. 

• The Fin Heat Transfer Rate 

Figures 5e and 6e show the comparison of the sur- 
face heat transfer rate profile against the experimental 
data on the fin surface at s = 0.167 and s = 0.1962 
m as a function of distance from the flat plate, for 
the 10° and the 15° fin angle, respectively. All the 
results with the two-equation models assuming fully 
turbulent flow from the fin leading edge overpredict 
the heat transfer rate by almost 175% in the laminar 
boundary layer. Close agreement is obtained when 
laminar flow is assumed on the fin wall. The transition 
Jc — u> model 21 overpredicts the heat transfer rate by 
almost 85%. The heat transfer rate on the fin surface 
near the corner inside the plate boundary layer is not 
well predicted. Further studies on transition and 3-D 
combination of laminar and turbulent boundary layers 
is needed. A simulation assuming turbulent boundary 
layer inside the plate boundary layer and transition 
modeling above the plate boundary layer edge did not 
improve the results. 

• Pitot Pressure Profiles 

Figures 5f and 6f show the comparison of the pitot 
pressure profiles against the experimental data as a 
function of distance from the flat plate. For the 10° 
fin angle, the data is at x = 0.1723 m and z — 0.064, 
0.0127, 0.0191, 0.0254, 0.0381, and 0.0508 m. The re- 
sults obtained with the standard k - w model show 
close agreement with the experimental data. Small 
differences are observed near the intersection of the 
boundary layer edge and the oblique shock wave. For 
the 15° fin angle, the data is at x = 0.1744 m and z = 
0.0191, 0.0254, 0.0381, 0.0508, 0.0635, and 0.0762 m. 
The results obtained with the baseline k — w model 
show close agreement with the experimental data, ex- 
cept near the intersection of the boundary layer edge 
and the oblique shock wave. 

• Pitot Pressure Contours 

Figures 5g and 6g show the pitot pressure con- 
tours in the crossed section corresponding to the pre- 
vious Figures 5f and 6f, respectively. The contour lines 
show the strong oblique shock and the wave structure 
of the of the main crossflow yortex. A comparison 
of the pitot pressure contours against the experimen- 
tal data is shown in Figures 5h and 6h, respectively. 
The agreement is very good for the 10° fin angle case, 
although some differences are observed near the inter- 
section of the turbulent boundary layer edge and the 
oblique shock wave. Similar results are observed for 
the 15° fin angle case, except for slightly larger differ- 
ences in the crossflow vortex. 

CONCLUDING REMARKS 

A 3-D Navier-Stokes numerical simulation of the 
3-D shock wave/ boundary layer interaction with a 
10° and a 15° fin angle has been conducted. Compar- 
isons between the results obtained with the baseline 
jfc — cj, the standard k — w, and the Baldwin-Lomax 
algebraic model agree in general with the experimen- 
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tal data. The simulations with the baseline Jfc — u 
turbulence model provides excellent agreement with 
the adiabatic flat plate data, good agreement with the 
cooled flat plate data, and fairly good agreement with 
the complex shock wave/ boundary layer interaction 
data. The flow structures are fairly well capture, in 
particular, the external planar oblique shock and the 
shock-wave/ boundary layer interaction zone are free 
of oscillations. There is a need to improve the models 
in order to avoid the underprediction of heat transfer 
rate in the crossflow separation zone. The model of 
“conical” approximation of the flow is not supported 
by the present simulations. Finally, the present numer- 
ical simulations are very efficient and require less cpu 
time than previous simulations with other algorithms, 
which is an important requirement to test different 
turbulence models in complex flows. 
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Figure la. 3x61x61 flat plate grid. 
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Figure 2a. 91x61x3 flat plate grid. 
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Figure lb. Skin friction on insulated plate. 
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Figure 2b. Skin friction on cooled plate. 
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Figure lc. Velocity profile on insulated plate. 
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Figure 2c. Velocity profile on cooled plate. 
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Figure Id. Velocity profile on insulated plate. 
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Figure 2d. Velocity profile on cooled plate. 
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Figure 3. 61x41x61 SWBLI grid showing plate, fin ,and outflow planes. 



Figure 4c. Inflow SWBLI temperature profile. Figure 4d. Inflow SWBLI density profile. 
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Figure 5b. Comparison of pressure profile on 
fin surface at s=0.1834 m from fin leading edge. 



Figure 6b. Comparison of pressure profile on 
fin surface at s=0.1834 from fin leading edge. 
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Figure 5c. Comparison of skin friction profile 
on plate surface at x=0.155 m. 
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Figure 6c. Comparison of skin friction profile 
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Figure 5d. Comparison of heat transfer rate Figure 6d. Comparison of heat transfer rate 

profile on plate surface at x=0. 1645 m. profile on plate surface at x=0.1645 m. 
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Figure 6e. Comparison of heat transfer rate 
-ome on fin surrace at s=0. 167 m and s=0. 1962 
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figure6f. Comparison of pitot pressure profi- 
les at x=0.1744 m. 



Figure 6g. Pitot pressure contours. 



Figure 6h. Comparison of pitot pressure 
contours at x=0.1/44 m. 
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